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I Abstract. Given a self adjoint operator A on a Hilbert space, suppose that that one 

QQ I wishes to compute the spectrum of A numerically. In practice, these problems often 

T-H I arise in such a way that the matrix of A relative to a natural basis is "sparse" . For 

example, doubly infinite tridiagonal matrices are usually associated with discretized 
I second order differential operators. In these cases it is easy and natural to compute 

■ the eigenvalues of large n x n submatrices of the infinite operator matrix, and to 

I hope that if n is large enough then the resulting distribution of eigenvalues will give 

. a good approximation to the spectrum of A. 

' While this hope is often realized in practice it often fails as well, and it can fail 

in spectacular ways. The sequence of eigenvalue distributions may not converge as 
n — > oo, or they may converge to something that has little to do with the original 
. operator A. At another level, even the meaning of 'convergence' has not been made 

' precise in general. In this paper we determine the proper general setting in which 

I , one can expect convergence, and we describe the asymptotic behavior of the n x n 

^ ' eigenvalue distrubutions in all but the most pathological cases. Under appropriate 

^ , hypotheses we establish a precise limit theorem which shows how the spectrum of A 

f ^ ' is recovered from the sequence of eigenvalues of the n x n compressions. 

• • , In broader terms, our results have led us to the conclusion that numerical prob- 

lems involving infinite dimensional operators require a reformulation in terms of 
C*-algebras. Indeed, it is only when the single operator A is viewed as an element 
' of an appropriate C*-algebra A that one can see the precise nature of the limit of 

the n X n eigenvalue distributions; the limit is associated with a tracial state on A. 
Normally, A is highly noncommutative, and in the main applications it is a simple 
C*-algebra having a unique tracial state. We obtain precise asymptotic results for 
the case where A is the discretized Hamiltonian of a one-dimensional quantum system 
with arbitrary potential. 
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1. Introduction. 

There are efficient algorithms available for computing the spectrum of self-adjoint 
n X n matrices. Using a good desktop computer, one can find all of the eigenvalues 
of a 100 X 100 Tocplitz matrix in a few seconds. It is natural to ask how these finite- 
dimensional methods might be applied to compute the spectrum of the discretized 
Hamiltonian of a one dimensional quantum system, or a self-adjoint operator in an 
irrational rotation C*-algebra. Considering the burgeoning capabilities of modern 
computers, it is unfortunate that there are apparently no guidelines available which 
indicate how one might proceed with such a project except in a few isolated cases. 
For example, we were unable to determine from the literature if such a program is 
even possible in general. 

The purpose of this paper is to develop some general methods for computing 
the spectrum (more precisely, the essential spectrum) of a self-adjoint operator in 
terms of the eigenvalues of a sequence of finite dimensional matrix approximations. 
We will see that there is a natural notion of degree of an operator (relative to a 
filtration of the underlying Hilbert space into finite dimensional subspaces) which 
generalizes the classical idea of band-limited matrices. More significantly, we will 
find that this kind of numerical analysis leads one to a broader context, in which 
the operator is seen as an element of a C*-algebra which is frequently simple and 
possesses at least one tracial state. The limit distributions arising from the finite 
dimensional eigenvalue distributions correspond to traces on this C*-algebra. 

Our intention is not to discuss algorithms per se, but rather to present an ef- 
fective program for utilizing the finite dimensional techniques to compute infinite 
dimensional spectra, along with a precise description of a broad class of operators 
for which it will succeed. Thus, we deal with issues arising from operator theory 
and operator algebras; we do not address questions relating directly to numerical 
analysis. 

In more concrete terms, given an orthonormal basis {ei,e2, . . .} for a Hilbert 
space H and an operator A e B{H), we may consider the associated infinite matrix 

{aij) 

ttij =< Aej, >, i, J = 1, 2, . . . , 

and one can attempt to approximate the spectrum of A by calculating the eigen- 
values of large n x n submatrices 



/ an 
\ani 



ai2 

0'n2 



■ ■ ain\ 

■ ■ <l2n 



This program fails dramatically in general. For example, if A is the simple unilateral 
shift (defined on such a basis by Ae^ = e^i+i, n > 1), then the spectrum of A is the 
closed unit disk and the essential spectrum of A is the unit circle, while oi^A^^ — {0} 
for every n because each A^ is a nilpotent matrix. 

Examples like this have caused operator theorists to view matrix approximations 
with suspicion, and to routinely look elsewhere for effective means of computing 
spectra. On the other hand, since the 1920s there has been evidence that for self- 
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Consider for example the foUowing theorem of Szego. Let / be a real function 
in L°°[— TT, tt] and let A be the multiplication operator 



Ai{x) = f{x)i{x), 



eeL^[-7r,7r]. 



With respect to the familiar orthonormal basis {e^ : n e Z}, en{x) = e**^^, the 
matrix of ^ is a doubly-infinite Laurent matrix whose principal n x n submatrices 
are finite self-adjoint Toeplitz matrices 



/ ao a_i 
ai ao 



a-n-l-2 



ao 



the numbers being the Fourier coefficients of /. Let {Ai, A2, . . . , A„} be the 
eigenvalues of repeated according to multiplicity. Szego's theorem asserts that 
for every continuous function tt : R — > M we have 



1.1 



lim -(«(Ai) + u{X2) + ■■■ + u{Xn)) = ^ [ u{f{x)) 

re— >oo n Ztt J _^ 



dx. 



See [8] [10]; a discussion of Szego's theorem can be found in [11], and [6] contains a 
meticulous discussion of Toeplitz operators. 

Let us reformulate Szego's theorem in terms of weak*-convergence of measures. 
Let a and h be the essential inf and essential sup of / 

a — ess inf f[x) 

— ■K<X<TT 

b = ess sup f{x), 

—n<x<n 

and let fi be the probability measure defined on [a, b] by 

1 



27T 



m{x : f{x) e S}, 



m denoting Lebesgue measure on [— 7r,7r]. Then 1.1 asserts that the sequence of 
discrete probability measures 

P'n = -((^Ai +(^A2 + H5a„), 

Sx denoting the unit point mass at A, converges to n in the weak*-topology of the 
dual of C[a,b], 



lim 

n—KX> 



u{x) djJinix) 



u{x)diJ,{x), uEC[a,b]. 



The spectrum of A coincides with its essential spectrum in this case, and is exactly 
the closed support of the measure /x. 

Once one knows that the measures converge in this way to a measure sup- 
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at which the eigenvalues of accumulate at points of the essential spectrum. For 
example, suppose that we choose a point A e cre(^). Then for every open interval 
/ containing A wc have fi{I) > 0. If wc choose positive numbers a and /? which are 
close to and satisfy a < < (3, then in the generic case where A has no 
point spectrum we may conclude that 

na < ij^{a{An) d I) < n/S 

when n is sufficiently large. Of course, the symbol ^ means the number of eigenval- 
ues in the indicated set, counting repeated eigenvalues according to the multiplicity 
of their occurrence. In a similar way, the density of eigenvalues of A^ tends to be- 
come increasingly sparse around points of the complement of the essential spectrum. 

In this paper we will show that, while convergence of this type cannot be expected 
in general, it does persist in a much broader context than the above. Actually, this 
work was initiated in order to develop methods for computing the essential spectrum 
of certain tridiagonal operators which are associated with the discretized Hamil- 
tonians of one-dimensional quantum mechanical systems (see section 5). These 
operators serve to illustrate the flavor of our more general results. Starting with 
a real valued continuous function V of a real variable, we fix a real number 9 and 
consider the tridiagonal operator T defined in terms of a bilateral orthonormal basis 
{en-.ne Z} by 

TCfi = 6jT,_i -|- dfiCfi -\- Cfi+ii 

where the diagonal sequence is defined by dn = F(sin(n0)), n e Z. The compression 
Tn of T to the span of {ck : —n < A; < n} is a {2n+l) x {2n+l) tridiagonal matrix 
whose eigenvalues can be computed very rapidly using modern algorithms. For 
example, we show that if Hn is the discrete probability measure defined by the 
eigenvalue distribution of 

then there is a probability measure having closed support ae{T), such that for 
every continuous function / : R — > IR we have 



lim / f{x)diin{x) = / f{x)dii{x). 



See formula 5.1. 

In these cases the limiting measure /x can be identified in rather concrete terms. 
Briefly, if O/tt is irrational then the pair of unitary operators 

Uen = Cn+l 

Vp = e^'^^e 

generate an irrational rotation C*-algebra Ag. Ae has an essentially unique repre- 
sentation in which its weak closure is the IIi factor R. T becomes a self-adjoint 
operator in i?, and the measure ji is defined by 

..to\ I TT' row o 1^ la 
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where r is the normahzed trace on R and Et is the spectral measure of T. 

In sections 2 and 3 we consider a rather general question. Given a self-adjoint 
operator A and an orthonormal basis for the underlying Hilbert space, let (oy) be 
the matrix of A relative to this basis. How must one choose the basis so that the 
essential spectrum of A can be computed in principle from the eigenvalues of the 
finite dimensional principal submatrices of {aij)l Roughly speaking, we show that 
this is possible if the diagonals of the matrix (a^) tend to zero sufficiently fast. It 
is not possible in general. The family of all operators whose matrices satisfy this 
growth condition is a Banach *-algebra but it is not a C*-algebra. 

In section 4 we seek more precise results which relate to weak*-convergence 
as above. It is appropriate to formulate these results in terms of C*-algebras of 
operators. We introduce the notion of degree of an operator (relative to a filtration 
of the underlying Hilbert space into finite dimensional subspaces), and show that 
one can expect convergence of this type in the state space of the given C*-algebra 
when the C*-algebra contains a dense set of finite degree operators and has a unique 
tracial state. These results are applied in section 5 to the discretized Hamiltonians 
described above. 

2. Spectral asymptotics of sequences of matrices. 

Let ^ be a bounded self-adjoint operator on a Hilbert space H and suppose 
that we have a sequence of finite-dimensional subspaces of H with the property 
that the corresponding sequence of projections ~ converges strongly to 
the identity. Normally, the subspaces will be increasing with n, the main 
examples arising in the cases where we start with an orthonormal basis {e„ : n = 
1,2,...} (resp. {cji : n e Z}) and choose = [ei, 62, . . . , e„] (resp. Hn = 

For every n > 1, let An denote the compression 

2.1 An = PnA\H^. 

More generally, given an arbitrary sequence Ai, A2, . . . of finite dimensional self- 
adjoint operators, we introduce two "spectral" invariants A, Ag of the sequence 
as follows. A is defined as the set of all A G M with the property that there is a 
sequence Ai, A2, . . . such that A^ e (T{An) and 

lim An = A. 

n — »oo 

Notice that A is closed. Indeed, a real number A belongs to the complement of A 
iff there is an open set U containing A and an infinite sequence of integers ni < 
77-2 <.. . with the property that (t(A„j,) fl t/ = for every k = 1, 2, ... . Every point 
of U shares this property with A, and hence the complement of A is open. If the 
sequence {A^} is bounded in norm then A is compact. Of course, A is considerably 
smaller than the set 

L = HniaiAn) U a{An+i) U . . .) 

of limit points of the sequence a{Ai), a{A2), . . . , since limits along subsequences 
may not qualify for membership in A. 

For every n > 1 and every set S of real numbers, let Nn{S) be the number 
of eigenvalues of A^ which belong to S, counting multiple eigenvalues accord- 
ing to their multiplicity. A^n(') is an integer-valued measure which takes values 



6 



WILLIAM ARVESON 



Definition 2.2. 

(1) A point A e IR is called essential if, for every open set U containing A, we 
have 

lim Nn{U) = oo. 

n— >oo 

The set of essential points is denoted Ag. 

(2) A e ]R is called transient if there is an open set U containing A such that 

SUpNn{U) < OO. 
n>l 

Remarks. It is clear that Ag C A. Moreover, notice that A e M is nonessential iff 
there is an open set U containing A and an infinite sequence of integers ni < n2 < 
. . . such that 

Nn,{U)<M <oo 

for every k = 1,2,.... Thus the nonessential points are an open set in R, and we 
conclude that the set Ag of essential points is closed. 

It is conceivable that A may contain anomalous points which are neither tran- 
sient nor essential. In small neighborhoods U of such a point one would find sub- 
sequences nfc for which Nn^.{U) remains bounded, and other subsequences for 
which Nmk{U) is unbounded. Fortunately, in most situations one has a more man- 
ageable dichotomy (see Theorem 3.8). At this level of generality, one cannot say 
much. We do have the following 

Theorem 2.3. Assuming that the sequence Ai, A2, . . . arises from an operator A 
as in 1.1, then we have cr[A) C A and a^^A) C Ag. 

proof of a {A) C A. 

Suppose that A is a real number which does not belong to A. We will show that 
^ — Al is invertible. Since A ^ A, there is an e > and a subsequence ni < n2 < . . . 
such that 

(7(^„Jn(A-e,A + e) = 
for every A; > 1. Thus the distance from A to a{An^) is at least e, hence 

II (A,, -Al)-^ ||<l/e 

for every k. Let 

Bk = {An,-Xl)-^Pn,. 

Since Bi, S2, . . . is a bounded sequence of self-adjoint operators and the ball of 
B{H) of radius 1/e is weakly sequentially compact, {-Bfc} has a weakly convergent 
subsequence. Thus, by replacing the original sequence ni,n2,... with a subse- 
quence of itself, we may assume that there is a bounded operator B such that 

lim < Sfe^, 7/ >=< S^, ?7 >, C,V^H. 

k—>-oo 



Since 
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for every k and since — Al)P„j. converges to ^ — Al in the strong operator 
topology, it follows that the product {An,, — Al)S„j. converges weakly to {A — X1)B. 
Hence from 2.4 we may conclude that 

{A - X1)B = weak limP„^ = 1, 

fe— »oo 

proving that A — Al is right-invertible. Left-invertiblity follows by taking adjoints 
in the preceding equation. 

proof ofa^iA) C A^. 

Suppose that A ^ Ag. We will show that A — Al is invertible modulo trace 
class operators, and hence A ^ o^e{A). By the hypothesis on A, there is a sequence 
ni < 77-2 < . . . of positive integers and a pair of positive numbers e, M such that 

Nn^{X-e,X + e) < M 

for every k = 1,2,.... For each k, consider An^ to be a self-adjoint operator 
in B{Pnf,H) and let Ek be the spectral measure of this operator. Then Qk = 
Ek{X — e, X + e) is a projection commuting with A^,. satisfying dimQ^ < M, such 
that the restriction Bj. of An,. —XI to the range of — Qk is an invertible operator 
with II B^^ II < 1/e. 

By using compactness again and passing to a subsequence as in the preceding 
argument, we can assume that both sequences {B^^^P^^. —Qk)} and {Qk} converge 
in the weak operator topology 

B^\Pn,-Qk)^C 

Qk Q- 

Notice that Q is a positive trace class operator. Indeed, since the set 

{T e B{H) : T > 0, trace(T) < M} 

is closed in the weak operator topology, the assertion is evident. Finally, using the 
fact that An^. — XI converges in the strong operator topology to ^ — Al and 

(An, - Xl)B-\Pr,, - Qk) = Pn, - Qk 

for every k, we may take weak limits in the latter formula to obtain 

{A - X1)C = 1-Q, 

as required □ 

Remarks. The appendix contains an example which shows that the inclusions of 
Theorem 2.3 can be proper. 

One can imagine computing the eigenvalues of A^ for every n = 1,2,..., picking 
a point A G A, and observing the distribution of eigenvalues in the vicinity of A as n 
becomes large. If A is essential then the number of eigenvalues in any small interval 
containing A will increase without limit. If A is transient, then there are positive 
integers p < q with the property that eventually you always see at least p points 
in the interval but never more than q. Moreover, after p and q are appropriately 
adjusted, then the same behavior will occur in every sufficiently small interval 
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3. Filtrations and their Banach algebras. 

Let ^ be a self adjoint operator on a Hilbert space H. Suppose we have an 
orthonormal basis for H, indexed either by N or Z, and we form the matrix (aij) 
of A relative to this basis. If we form a sequence of finite-dimensional compressions 
of (ay) then we may compute A and Ae in principle, and the general results of 
section 2 imply that (Jf^{A) C Ag. In this section we find conditions on the matrix 
(oy) which guarantee first, that every point of A is either essential or transient 
and second, that a"e(A) = Ag. The criterion is that the series 

Efc|fc|^/^4 should 

converge, where is the sup norm of the kih. diagonal of {aij). We conclude this 
section with some comments about the nature of transient points. 

In order to deal with both unilateral and bilateral orthonormal bases as well as 
more general situations in which the dimensions of compressions increase in uneven 
jumps, it is necessary to work with filtrations. We introduce the concept of degree of 
an operator and a Banach *-algebra D{T) of operators associated with a filtration 

It is important for the applications that it should be easy to estimate the norm 
in D{J^) (see Proposition 3.4). 

Definition 3.1. 

(1) A filtration of if is a sequence = {Hi, H2, . . .} of finite dimensional 
subspaces of H such that C H^^i and 



(2) Let JF = {Hn} be a filtration of H and let P„ be the projection onto H^. 
The degree of an operator A e B{H) is defined by 

deg(^) = suprank(P„^ - APn). 

n>l 

Remarks, deg(^) is either a nonnegative integer or +00. Moreover, we have 

deg(A*) = deg(A), deg(AA) = deg(A) for all A ^ 0, 

deg{A + B) < deg{A) + dcg(B), dcg(AS) < deg(A) + deg(S). 

Only the last of the four properties is not quite obvious; it follows from the fact that 
the map A~^[Pn,A]^ PnA-APn is a derivation, [P^, AB] = [P^, A]B + A[Pn,B], 
which implies 

rank[Pn, AB] < rank[Pn, A] + rank[Pn, B]. 

We conclude from these observations that the set of all finite-degree operators in 
B{H) is a self-adjoint unital subalgebra of B{H). 

Remark. Notice that for any operator A and any projection P we have 

PA-AP = PA{1 - P) - (1 - P)AP, 

hence 

rank(PA - AP) = rank(PA(l - P)) + rank((l - P)AP) 

r>\ A*r>\ i i./'/'-i r>\ A r>\ 
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It follows that deg(A) < oo iff both operators B = A and B = A* satisfy the 
condition 

suprank((l — Pn)BPn) < +00. 

n>l 

Thus, operators of finite degree are abstractions of band-limited matrices. 

We associate a Banach *-algebra D{J^) to a filtration in the following way. 
Let D{T) denote the set of all operators A e B{H) having a decomposition into an 
infinite sum of finite degree operators Af- 

00 

(3.2) A = Y,Ak 

in such a way that the sum 

00 

s = 5^(l + deg(Afe)i/2)||Afc|| 

/s=l 

is finite. Notice that finiteness of s ensures that the series 3.2 converges absolutely 
with respect to the operator norm. We define \A\j: to be the infimum of all such 
sums s which arise from representations of A as in 3.2. 

Proposition 3.3. With respect to the norm \ ■ \jr and the operator adjoint, D[T) 
is a unital Banach *-algebra in which \l\jr = 1. 

proof. Let V be the *-algebra of all finite degree operators in B{H), and consider 
the function : D — > R"'' defined by 

<P{A) = {l + idegAy/^)\\A\\. 

The norm on D{J^) is defined by 

|^|^ = inf^0(ylfc), 

k=l 

the infimum extended over all sequences Ai, A2, . . . inV satisfying 

00 

k=i 

together with '^j^<t>{Ak) < +00. The fact that {D{J^), \ ■ is a Banach *-algebra 
follows from straightforward applications of the following properties of 0: 

(1) HA) > \\A\\. 

(2) (j){XA) = |A| ■ Pll, A e C. 

(3) (^(A*)=</>(A). 

(4) cPiAB) < <P{A)<p{B). 

For example, (4) follows from the fact that the weight sequence wo,wi,W2, . . . 
defined hy Wk — 1 + k^^'^ satisfies 

Wk+j < WkWj, k,j > 0, 

together with 

deg(^5) < deg^ + deg5. 

The identity of B{H) is of degree 0, hence \l\jr < 1, while the inequality |1|^ > 1 
is valid for the unit of any Banach algebra □ 

D{J^) is certainly not a C*-algebra, but it is dense in B{H) in the strong operator 
topology. The following gives a concrete description of operators in D{T) for one 
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Proposition 3.4. Let {e„ : n e Z} he a bilateral orthonormal basis for a Hilbert 
space H and let T = {Hi, H2, . . .} be the filtration defined by 

Let {ttij) be the matrix of an operator A G B{H) relative to {en}, and for every 
keZ let 

dk = sup \ai+k,i\ 

be the sup norm of the kth diagonal of (aij). Then 

+00 

(3.5) 1^1^ < J2 (1 + |2A;r/')c/fc. 

k= — oo 

In particular, A will belong to D{T) whenever the series \k\^/'^dk converges. 

proof. We may assume that the sum on the right side of 3.5 converges. Let Dfe 
be the operator whose matrix agrees with (ajj) along the kth diagonal and is zero 
elsewhere. Notice that ||-Dfc|| = dfe, and hence 

+00 
fe=— 00 

the sum on the right converging absolutely in the operator norm because 

Y,\k\^'^dk<oo. 

k 

We claim that degD^ < 2\k\ for every k E To see that, fix n = 1,2,... and 
let Pn be the projection onto H^. The operator maps the range of Pn into the 
range of Pn+\k\, and a careful inspection of the two cases k > and k < shows 
that the dimension of (1 — Pn)Pn+\k\ is at most |A;| for every n > 1. It follows 
that the rank of either operator (1 — Pn)DkPn or (1 — Pn)DlPn is at most 
independently of n. Thus from a previous remark we conclude that degD^ < 2\k\. 

Now ll-Dfcll = dk, and since the series |A;|-^/^cifc converges we must have 
X)fc ll-Dfcll < 00. Hence 

+00 

A= J2 

k= — oo 

is an absolutely convergent sum of operators satisfying degD>fc < 2|/c|, from 
which 3.5 follows □ 

In particular, any operator whose matrix (0^) is hand-limited in the sense that 
aij = whenever \i — j\ is sufficiently large, must belong to D{T). 

The following estimate occupies a key position, and it explains why the norm on 
-D(J^) is defined as it is. For B e B{H), we write ||S||2 = (trace(S*S))V2 for the 
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Lemma 3.6. Let {Hi,H2, . . .} be a filtration of H and let P„ he the projection 
onto Hn- Then for every A e D{J^) we have 

sup \\APn - PnAPnh < \A\jr < OO. 
n>l 

proof. We claim first that for every A e B{H) and every n = 1, 2, . . . we have 
(3.7) trace] AP„ - Pn^i'nP < deg{A)\\Af, 

where for an operator B e B{H), \B\ denotes the positive square root of B*B. 
Indeed, since tracejPp < rank(P)||P|p, we can estimate the left side of 3.7 as 
follows 

trace|(l - Pn)^Pn|^ < rank((l - P^)AP^))\\{1 - Pn)APnf 

= rank((l - Pn){AP^ - PnA))||(l - Pn)APnf 
< vankiAPn - PnA)\\Af < deg{A)\\Af. 

By taking the square root of both sides of 3.7 we obtain 

\\APn-P,,APj2<deg{Ay/'\\A\\. 

Now suppose A e D{T), and choose finite degree operators Ak such that 

5](l + deg(^fe)V2)p,||<oo 

k 

and A = Ai-\-A2 + Using the triangle inequality for the Hilbert-Schmidt norm 

we see from 3.7 that 

CXD 

\\PnA- PnAPnh < 5^deg(Afc)i/2||ylfc|| 

OO 

< J](l + deg(^fc)V2)p,||. 

/s=l 

The desired inequality follows by taking the supremum over n and the infimum 
over all such sequences ^i, ^2, . . . □ 

Remark. Notice that Lemma 3.6 implies that for every A e D{T) we have the 
following uniform estimate of commutator 2-norms 

sup \\APn - PnA\\2 < V2\A\jr < +00. 
n 

Now let A be a self-adjoint operator in B{H), let = {Hi, H2, . . . } be a filtra- 
tion, and let A^ be the compression of A to Hn- The following result is our basic 
description of the spectrum of A in terms of the eigenvalues of It asserts 

that the essential spectrum of A is precisely the set of essential points, and that all 
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Theorem 3.8. Assume that A = A* belongs to the Banach algebra D{T). Then 

(i) a^{A) = Ae. 

(ii) Every point of A is either transient or essential. 

proof. We will prove that every point in the complement of the essential spectrum 
of A is transient. Notice that both statements (i) and (ii) follow from this. Indeed, 
because of Theorem 2.3 we know that ae{A) C Ag. On the other hand, since 
no transient point can be an essential point, the above assertion implies that the 
complement of ae{A) is contained in the complement of Ae; hence (i). Armed with 
(i), we see that A \ Ag is the same as A \ ae{A), which by the assertion is contained 
in the transient points of A. Thus we obtain (ii). 

To prove this assertion, choose a point A in the complement of ae{A). In order to 
show that A is transient, we claim first that there is a finite-dimensional projection 
Q and a positive number e such that for every real number e [A — e, A + e] , the 
operator 

A + Q- nl 

is invertible. Indeed, since A ^ cre(A) it follows that either A is an isolated point of 
the spectrum of A whose eigenspace is finite dimensional, or else ^ — Al is invertible. 
In either case 

{^eH:A^ = AC} 

is finite dimensional and if Q denotes the projection onto this eigenspace, then Q 
commutes with A and ^ + Q — Al is an invertible self-adjoint operator. Because 
the invertible operators are an open set it follows that there is a neighborhood U 
of A with the property that A + Q — /j,! is invertible for every /j, E U. We can take 
[A — e, A-h e] to be an appropriate closed interval about A which is contained in such 
a neighborhood. 

Note too that by continuity of inversion, there is a positive number M such that 

(3.9) II (A + Q-//l)-^ ||< M 

for every satisfying ||U — A| < e. 

We will show that if V is the open interval (A — e, A + e) and N^iV) is the number 
of eigenvalues of An which lie in V, then the sequence Nn{V) is bounded. Thus A 
must be a transient point, and that will complete the proof. To that end, note first 
that by Lemma 3.6 there is a positive number N (in fact, N = |^|^ will do) such 
that 

trace|APn - PnAPnf < N, 

for every n = 1, 2, . . . . Fix n, and let Ai, A2, . . . , Ap be the set of all eigenvalues of A^ 
which lie in V, with repetitions according to the multiplicity of the corresponding 
eigenspace. Actually, we will show that 

(3.10) p<M^N + diinQ. 

Since the right side of 3.10 does not depend on n it will follow that A is transient, 
completing the proof. 

To prove 3.10, we choose for each k = 1,2, ... ,p a, unit vector ^fe in the domain 
of An (i.e., the range of Pn) such that An^k = ^k^k, and such that {^1, . . . ,^p} 
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different eigenvalues of a self-adjoint operator are mutually orthogonal; in the case 
of a multiple eigenvalue we simply choose an orthonormal basis for that eigenspace. 
Note first that for each k between 1 and p, 

(3.11) A^k - AfcCfe = (1 - Pn)A^k- 
Indeed, since Pn{A^k - ^k^k) = ^n^fc - ^k^k = 0, we can write 

A^k - AfcCfc = (1 - Pn){A^k - Afc^) = (1 - Pn)A^k- 

We claim now that for each k between 1 and p, we have the inequality 

(3.12) 1 < M2 II (1 - P^)A^k ||2 + < Q^^, > . 

To see that, write 

1- < Qa,e/c > = II (1 - Q)a II' < II + Q - Afc)(l - Q)a IP 

= II (A - Afcl)(l - Q)Ck f = II (1 - QXA^k - AfcCfe) II 
= II (1 - Q)(l - Pr,)ACk f < II (1 - P^)Aa f . 

The first inequality follows from 3.9, the identity {A — Xk) {1 — Q)^k = {^ — Q){A^k — 
Afc^fc) follows from the fact that Q and A commute, and the last equality and 
inequality follow respectively from 3.11 and the fact that 1 — Q is a contraction. 
Summing the inequality 3.12 on k we obtain 

p V 

P < M2 ^ II (1 - Pr.)Aik f + J] < Qik. ik > 

k=l k=l 

< trace|(l - P^)^P^|^ + trace(Q) 

< M'^N + dimQ, 



and 3.10 is established □ 

Remark 3.13. We conclude this section with some remarks about transient points 
and other spurious eigenvalues. To provide a context for these remarks, let us 
return to the classical setting of Szego's theorem, in which Ai,A2,... is a sequence 
of Toeplitz matrices. Let / be a bounded measurable real- valued function on the 
interval [— 7r,7r], and consider the Fourier series of / 
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For n = 1, 2, . . . let A^ be the Toeplitz matrix 
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Let a and h be respectively the essential inf and sup of / over [tt, tt] . It is known that 
the eigenvalues of the sequence {An : n = 1,2,...} tend to fill out the entire interval 
[o, b] (sec [11, pp 201-202] and [12, Lemma ILl]). In more concrete terms, suppose 
that A G (a, b) and that we are looking in a small neighborhood U = (A — e, A + e) 
of A. Then U fl a (An) 7^ for sufficiently large n, regardless of where A is located 
in the interval (a, b). The point we want to make is that this can be misleading, 
and that what one is observing here may have nothing to do with the spectrum of 
the operator A. 

Indeed, for these operators A we have 

a{A) = ae{A) = R, 

where R is the essential range of the function /. Suppose that / is chosen so that 
its essential range is disconnected. If we choose a point AG [a, 6] \ i? then for every 
e > we will find the sets of eigenvalues (A — e, A + e) fl cr(A„) to be permanently 
nonempty for large n, in spite of the fact that (A — e, A + e) contains no points of 
the spectrum of A if e is small. 

Returning to the context of Theorem 3.8, it is not hard to give examples of band- 
limited matrices (0^) which exhibit the same kind of behavior, namely (j{A) = 
ae{A) while A contains points not in the spectrum of the operator. In this case, 3.8 
implies that such points will be transient. More generally, if one starts with a band- 
limited operator whose spectrum is properly larger than its essential spectrum, then 
every point of <j{A) \ ae{A) will be a transient point (Theorem 3.8). Of course, for 
such an operator there may also be other transient points which do not belong to 
(7(^). Unfortunately, we do not know how to distinguish between transient points 
which belong to the spectrum and transient points which do not. Perhaps it is not 
even possible to do so. In any case, these observations have led us to the conclusion 
that transient points should be ignored. In doing that one is also ignoring points 
of o'{A) \ ae{A); but since such points are merely isolated eigenvalues having finite 
multiplicity, the cost is small. 

We conclude that if one is interested in approximating spectra of infinite dimen- 
sional operators using matrix truncations, then one should not simply look for the 
eigenvalues of the truncations. Rather, one should weight the count of eigenvalues 
in a way which eliminates spurious ones occuring in the vicinity of transient points. 
The simplest example of such weighting is the sequence of 'densities' 

1 

n 

where N^iU) is the number of eigenvalues of the nxn truncation of A which belong 
to the open set U. Indeed, in the present context we obviously have 

lim^ = 

n—>oo Tl 

whenever U is a, sufficiently small neighborhood of a transient point A, simply 
because the sequence Ni{U), N2{U), ... is bounded. If A is an essential point on 
the other hand, then we require conditions under which the limit 
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exists and is positive for every open neighborhood U of A. The existence of this 
limit in the case of the multiphcation operators A discussed above is a consequence 
of Szego's theorem, as we have pointed out in the introduction. In the following 
section we take up the asymptotic behavior of these densities in a more general 
context. 

4. Filtrations and convergence of densities. 

In this section we shift our point of view somewhat. Rather than consider single 
self-adjoint operators we consider a concretely presented C*-algebra of operators 
A C B{H) and we introduce the concept of an ^-filtration. This is an abstraction 
of an orthonormal basis for H with respect to which a dense set of the operators in 
A are band-limited. The point we want to make is that in order to use the methods 
of this paper to compute the spectrum of a self-adjoint operator in a C*-algebra C, 
one should first look for a faithful representation tt : C — > B{H) with the property 
that the C*-algebra A = 7r(C) admits a natural ^-filtration, and then proceed with 
the analysis associated with Theorem 3.8 or Theorem 4.5 below. 

After presenting some examples, we prove a counterpart of Szego's theorem which 
is appropriate for C*-algebras. In order to keep the statement and proof as simple 
as possible, we assume that the C*-algebra has a unique trace. While this is not 
the most general formulation possible, it does provide the basis for the applications 
which will be taken up in the next section. It will be convenient (though not 
necessary) to assume that C*-algebras of operators on H contain the identity of 
B{H). 

Definition 4.1. Let A C B{H) be a C*-algebra. An ^-filtration is a filtration 
of H with the property that the *-subalgebra of all finite degree operators in A is 
norm-dense in A. 

Example 1. The simplest example is that in which A is an AF-algebra having a 
cyclic vector ^. If .Ai C ^2 ^ . . . is a sequence of finite dimensional *-subalgebras 
of A whose union is norm-dense in A, then Hn = [AnC] defines a filtration for which 
every element of the union UnAn has finite degree. 

Example 2. Suppose on the other hand that we have an action of Z on a compact 
Hausdorff space X for which there is a point xq E X whose orbit under Z is dense. 
Then one can write down a representation tt of the C*-algebraic crossed product 

C = Z X C{X) 

and a natural filtration for the C*-algebra A = 7r(C). To see that, let {cn : n e Z} be 
an orthonormal basis for a Hilbert space H. The representation ttq : C{X) — > B{H) 
and the unitary operator U e B{H) defined by 

7i'o(/)en = f{n ■ XQ)en 

Uen= Cn+l 

define a covariant pair for the induced action of Z on C{X), and hence the pair 
(tto, U) defines a representation tt of the crossed product. Let 
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and let V be the image of C{X) under ttq. It is clear that every element of V has 
degree zero, and that deg(t/) = 1. Thus any operator expressible as a finite sum 

n 

k=—n 

where G V, has finite degree. Since these operators are dense in A it follows 
that {Hi C ^ • • • } is an ^-filtration. 

Example 3. Let T be the Toeplitz C*-algebra. We represent T as the C*-algebra 
generated by a simple unilateral shift S, acting on an orthonormal basis {ei, 62, . . . } 
for H via — Cn+i- Let J' — {H^i : ??. > 1} be the filtration H^i — [ei, 62, . . . , Cn]- 
Notice that deg(5') = 1. Indeed, letting Pn be the projection onto Hn, we have 

^Pfi — PnSPn + {Pn+1 ~ Pn)SPm 

while PnS = PnSPn- Hencc SPn — PnS = (Pn+1 — Pn)SPn is a rank-one operator 
for every n. It follows that the *-algebra Tq generated by S consists of finite degree 
operators, and Tq is dense in T. Thus .7-" is a T-filtration. 

We require the following estimate. 

Proposition 4.2. Let T = {iin\ be a filtration of H, let Pn be the projection onto 
Hn, and let Ai, A2, . . . ,Ap be a finite set of operators in B{H). Then for every 
n — 1,2, . . . we have 

p 

trace|Pn^l^2 • --ApPn - PnAiPnA2Pn ■ . . PnApPn\ < \\Ai\\ . . . || ^ degAfc. 

k=l 

proof. Without loss of generality, we can assume that each operator Ak is a con- 
traction and has finite degree. Let us fix n, and define 

Ap = P^AiA2 ■ ■ ■ ApPn — PnAiPn,A2Pn ■ ■ ■ PnApP-n- 

We claim first that there are contractions Bi, B2, ■ ■ ■ , Bp,Ci,C2, ■ ■ ■ ,Cp such that 

p 

Ap = J2Pkil - Pn)AkPnCk. 
k=l 

This is easily seen by induction on p. It is trivial for p = 1 (take B = C = 0). 
Assuming it valid for p, we have 

^p+l = -Pl^l • • -^^(1 — Pn)Ap+iPn + ApAp+iPn, 

and the conclusion is evident from the induction hypothesis on Ap. 
Each operator (1 — Pn)AkPn is of finite rank, since 

rank((l - Pn)ApPn) = rank((l - Pn){ApPn - P^Ap)) 
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Thus we may write 

p 

trace] Apl < ^trace|S/,(l — Pn)AkPnCk\ 

k=l 

p 

(4.3) <Y,\\Bk\\\\Ck\\tr^ce\{l - Pn)AkPn\ 

k=l 

P 

< ^trace|(l-P„)AfePn|, 
fe=i 

where we have used the elementary fact that for any finite rank operator F, 

trace|5FC| < ||S|| ||C|| trace|F|. 
Moreover, for such an F we have 

trace|F| < ||F|| • rank(F). 

Hence 

trace! (1 - Pn)AkPn\ < rank((l - Pn)AkPn) = rank((l - Pn){AkPn - PnAk)) 

< rank(^fcPn - PnAk) < deg^^. 

Utilizing the latter inequality in (4.3) we obtain 

p 

trace|Ap| < ^^degA/-, 

k=l 

as required. □ 

Let ^ be a concrete C*-algebra. A tracial state of ^ is a positive linear functional 
p satisfying p(l) = 1 and p{AB) = p{BA) for every A, B E A. One might try to 
construct tracial states of A by using a filtration — {Hn : n > 1} for H in the 
following way. Letting P„ be the projection on H^, we can define a state pn of 
B{H) by way of 

dim{Hn) 

The restriction of to PnB{H)Pn is a trace. Of course, pn does not restrict to 
a trace of A, but one might attempt to obtain a trace by taking weak*-limits of 
subsequences of {pn \a' ^ 1}; or their averages. In general, this program will 
fail. It will certainly fail if A is isomorphic to a Cuntz algebra On, n = 2,3, . . . , oo, 
since has no tracial states whatsoever. On the other hand, the following result 
shows that this procedure will succeed when JF is an ^-filtration. 

Proposition 4.4. Suppose that {Hi C H2 C . . . } is an A-filtration. For every 
n = 1,2, . . . put dn — (iim.Hn, and let pn he the state of A defined by 



qJA) = — trace(P„.4). 
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Let Rn be the weak* -closed convex hull of the set {p„, Pn+i, Pn+2, ■ ■ ■}■ Then HnRn 
is a nonempty set of tracial states. 

proof. i?i,i?2, • • • is a decreasing sequence of nonvoid compact convex subsets of 
the state space of A, and hence the intersection i?oo is a nonvoid compact convex 
set of states. We have to show that every element a G i?oo is a trace. Since the 
finite degree elements are dense in A and a is bounded, it suffices to show that 

a{AB) = a{BA) 

for operators A, B e A having finite degree. Because of the inequality 
\a{AB) - a(BA)\ < limsup \pn{AB) - pn{BA)\, 

n — >oo 

it suffices to show that \pn{AB) — pn{BA)\ ^ as n — > oo. But since 

iTax:e{PnAPnBPn) = tiace{PnBPnAPn) 
for every n, we may use the case p = 2 of Proposition 4.2 as follows 

\Pn{AB) - PniBA)\ = ^\{trace{PnABPn) - trace (P„SAP„)| 

dn 

< ^ {\tr 8iCe{PnABPn - PnAPnBPn)\ + \iTme{PnB AP^ - PnBPnAPn)\) 
dn 

<^\\A\\\\B\\idegA + degB). 

dn 

The right side obviously tends to as n ^ oo □ 

Remarks. Given a C*-algebra A Q B{H), it is reasonable to ask if it is always 
possible to find an ^-filtration. Proposition 4.4 implies that if an ^-filtration 
exists then A must have at least one tracial state. Thus, many concrete C*-algebras 
cannot be associated with a filtration. For example, if A is isomorphic to one of 
the Cuntz algebras On, n = 2, 3, . . . , oo, then A has no tracial states whatsoever, 
and hence ^-filtrations do not exist. 

We assume throughout the remainder of this section that A has a unique tracial 
state T. Then every self-adjoint operator A e A determines a natural probability 
measure p,A on R by way of 

/ + 00 
f{x)dpA{x)^T{f{A)) 
-oo 

for every / e Co(R). Clearly the closed support of pA is contained in the spectrum 
of A, and if r is a faithful trace then we have 

support (//^) = o^{A). 

Here is a more explicit description of pA- The GNS construction applied to r gives 
rise to a representation of A on another Hilbert space in which the image of A 
generates a finite factor M. The spectral measure Ea of A takes values in the 
projections of M. Letting tr denote the natural normalized trace on M we have 

PAiS) = tr{EA{S)) 

for every Borel set S QM.. The measure pA will be called the spectral distribution 
of A. 
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Theorem 4.5. Let ha be the spectral distribution of a self-adjoint operator A e 
A, and let [a,b] be the smallest closed interval containing o^{A). For each n, let 
dn = dim Hn and let Xi, X2, ■ ■ . , \d„ be the eigenvalues of = PnA \ Hn, repeated 
according to multiplicity. Then for every f e C[a, 6], 

lim ^(/(Ai)+/(A2) + --- + /(AdJ)= f f{x)d^iA{x). 

proof. Let t„ be the state of B{H) defined by 

Tn{T) = -^trace(PnT). 

Noting that t„ restricts to the normahzed trace on PnB{H)Pn, Theorem 4.5 wiU 
become evident once we have estabhshed the foUowing two assertions. 

(i) For every / e C[a,b], we have 

hm \TnifiA))-TnifiPnAPn))\ = 0. 

n— >oo 

(ii) For every B e A, 

hm Tn{B)^T{B), 
n— »oo 

where r is the tracial state of A. 
To prove (i), notice that since the sequence of hnear functionals 

/ e C[a, b] ^ Tn{f{A)) - Tn{f{PnAPn)) 

is uniformly bounded (an upper bound for their norms is 2), it is enough to prove 
(i) for polynomials /; and by linearity, we may further reduce to the case where 
/ is a monomial, f{x) — for some p = 1,2, ... . Finally, since the sequence of 
p-linear forms 

Bn{Ti,T2, . . ., Tp) = Tn{TiT2 . . . Tp) — Tn{PnTiPnT2Pn ■ ■ ■ PnTpPn) 

is uniformly bounded (again, ||-Bri|| < 2), it is enough to prove that 

lim \Tn{AP)-Tr,{{PnAPnY)\=^ 

n— >oo 

for operators A having finite degree. But in this case. Proposition 4.2 implies that 

\Tn{A^) - rn{{PnAPnY)\ < ^WW^^M. 

and the right side obviously tends to as n — > 00. 

Assertion (ii) follows immediately from Proposition 4.4. For if Rn denotes the 
sequence of convex sets of 4.5, then their intersection must consist of the singleton 
{r}; hence the sequence of states 

B e A^ tJB) = ^tracefP„S) 
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has a unique weak* cluster point r, hence it must actually converge to r in the 
weak* topology □ 

Remark 4.. 6. Let A — A* & A. If one drops the hypothesis that A has a unique 
trace, then the averages of the eigenvalues of PnA \ Hn may fail to converge. 
However, one can push the method of proof to obtain other useful conclusions in 
certain cases. For example, if we assume that there are sufficiently many traces 
in the sense that for every nonzero B E A there is a tracial state r such that 
t{B*B) > 0, then one can obtain the following information about the density of 
eigenvalues in the vicinity of a point of the spectrum of A. Let Nn{U) denote the 
number of eigenvalues of PnA \ which belong to the open set U, and let A be a 
point in the spectrum of A. Then for every neighborhood C/ of A there is a positive 
number e such that 

ti )i all sufficiently large n. We omit the details since this result is not relevant to 
our needs below. 

5. Applications. 

Most one-dimensional quantum mechanical systems are described by a Hamil- 
tonian which is a densely defined unbounded self-adjoint operator on L^(]R), having 
the form 

Hfix) = -lf"ix) + Vix)fix) 

V : M — s> M being a continuous function representing the potential. If one wishes 
to simulate the behavior of such a quantum system on a computer, one first has 
to discretize this differential operator in an appropriate way. Secondly, one has to 
develop effective methods for computing with the discretized Hamiltonian. For our 
purposes, we interpret the second goal to be that of computing the spectrum. 

In [1], [2], we discussed the problem of discretizing the above Hamiltonians in 
such a way as to preserve the uncertainty principle. We argued that the appropriate 
discretization is a bounded self-adjoint operator of the form 

(5.1) H^ = Ap^ + v{Q,) 

where and are the bounded self-adjoint operators given by 

P.fix) = ^ifix + a)-fix-a)), 

Qaf{x) = -sm{ax)f{x). 
cr 

Here, cr is a small (positive rational) number representing the numerical step size. 

In even the simplest cases, one cannot carry out explicit calculations of the 
spectrum of the discretized Hamiltonians 5.1 (see [3], [4], [5], [7], [9]). On the other 
hand, one can apply the methods of the prceding sections in a straightforward 
manner to obtain the spectrum of H^^ as the limiting case of finite dimensional 



C* -ALGEBRAS AND NUMERICAL LINEAR ALGEBRA 



21 



To see this, let U and V be the unitary operators 

Vf{x) = f{x + 2a). 

Then 5.1 can be written in the form H^^ = aA + f31, where a and /? are real numbers 
and A is a bounded self-adjoint operator of the form 

(5.2) A = V + V* + v( — (U -U*)). 

2i 

Here, : M — > M is an appropriately rescaled version of V. Thus, we are interested 
in the spectrum of A. Of course, A belongs to the C*-algebra C*{U, V) generated 
by U and V. Following the program outlined in section 4, we will first find a faithful 
representation of C* {U, V) in which there is a compatible filtration, and then we 
will look at the eigenvalues fo the approximating sequence of finite matrices. Since 
U and V obey the commutation relation 

(5.3) FC/ = e^^^VF 

and since o" is a rational number, it follows that C*{U, V) is an irrational rotation 
C*-algebra. Thus in order to define a faithful representation, it is enough to specify 
a pair of unitary operators Ui and Vi which satisfy the same commutation relation 
5.3. For that, consider a Hilbert space spanned by an orthonormal basis {e„ : n e 
Z} and define Ui and Vi by 

The = p-^^^^^p 
ViCn = en+l- 

The pair {Ui,Vi} is irreducible, it obeys the required commutation relations, and 
A becomes the tridiagonal operator Vi + V* + D, where D is the diagonal operator 

Dcn = v{— sm{2a^n))em n e Z. 

Let A = C*(C/i,Vi). Then the sequence of subspaces Hi, H2, . . . defined by 
Hn — [e-n, C-n+i, • • • , 6^] is an ^-filtration (see example 2 of section 4). Since con- 
cretely represented irrational rotation C*-algebras never contain nonzero compact 
operators, the spectrum of A coincides with the essential spectrum of A. Moreover, 
by Theorem 3.8, the spectrum of A is the set of essential points associated with 
the sequence of matrices {A^ n > 1}, An being the compression of A to Hn- 
The remaining points of A are all transient points. Finally, since irrational rotation 
C*-algebras are simple and have a unique trace, we may conclude from Theorem 
4.5 that for every / e Co(M) one has the following convergence of the densities of 
eigenvalues, 

(5.1) lim ^^^(/(Ai) + • • • + /(A2n+i) = r fix) dfiA{x), 

wehre ha is the spectral distribution of A and {Ai, . . . , A2n+i} is the eigenvalue list 
of An. 

We remark that this representation is particularly convenient for numerical 
eigenvalue computations using the QL or QR algorithms, because the matrices 
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6. Appendix. Failure of ae{A) — Ag. In this section we give an example 
of a self-adjoint operator ^ on a Hilbert space spanned by an orthonormal basis 
{ei, 62, ■ ■ ■}, such that the relation of A to the filtration = [ei, 62, . . . , Cm], 
n — 1,2, . . . is pathological. We show that if An is the compression of A to Hn, 
then Ae contains points not in the spectrum of A, and in particular Ae 7^ ae{A). 

Specifically, ^ is a reflection (i.e., a self-adjoint unitary operator) whose essential 
spectrum is {— whereas e Ag. In the example H = i^(N), e„ is the 
unit coordinate funtion at n, and A will be the operator in /^(N) induced by a 
permutation tt of N satisfying = id. Let In denote the interval {1,2, ...,n}, 
and let denote the number of elements of a set S. Suppose further that the 
permutation tt has the following property 

(A.l) lim #(7r(/„)\/n) = +oo. 

n— »oo 

For each n > 1, let An be the compression of the operator A^{k) = ^{n{k)) to 
Hn = l'^{In)- Since A maps to 67^(^)5 it follows that for every k in the set 
Sn — {k & In '■ 7r(/c) ^ In}, we have A^e^ = 0. Since the cardinality of Sn tends to 
infinity, it follows that is an eigenvalue of every A^, and that the multiplicity of 
this zero eigenvalue tends to infinity. Hence, must belong to Ag. 

Thus it suffices to exhibit an order two permutation having the property A.l. 
Let E and O denote respectively the set of even and odd integers in N. In order to 
define a permutation tt with tt^ = id, it is enough to specify a bijection / : i? — > O; 
once we have / then we can define tt by 



TT{k) 



f{k),keE 
f-\k),keO. 



We will define a bijection f : E ^ O having the following property with respect to 
the initial intervals In' 

(A.2) lim #{f{Enln)\In) = +oo. 

n— >oo 

The resulting permutation will have the desired property A.l. 

To construct /, split E into two disjoint subsets Ei U E2, where Ei — 4N is the 
set of multiples of 4 and E2 is the rest. Both Ei and E2 are infinite of course. First 
define / on Ei by 

(A.3) f{k) = k^ + l. 

Let Oi — f{Ei), and put O2 = 0\0i. Noting that O2 is infinite (everything in Oi 
is congruent to 1 mod 16, hence O2 contains all the other odd numbers), we can 
complete the definition of / by choosing it to be any bijection of E2 onto O2. 

To prove A.2, it suffices to show that #(/(£^i fl In) \ In) tends to infinity with 
n. Noting that 

EiDln^ {4:k:l<k< [n/4]}, 

we see that 
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Thus every element in the subset of Ei n In defined by 

P^ = {4k:^<k<l} 
gets mapped by / into the complement of In- Hence 

n In) \ In) > #Pn > ^(1 - n"^/^). 

This estimate implies that ^{f{E fl /„,) \ In) oo, and A. 2 follows. 

Remark. Notice that this misbehavior cannot be corrected by averaging the eigen- 
values of An as we did in section 4. Indeed, if Nn{U) denotes the number of 
eigenvalues of An which belong to an open set U, then the above estimate shows 
that for every neighborhood U of we have 

U„,inf^>l>0, 

n^oo dimlln 4 

in spite of the fact that does not belong to the spectrum of A. 
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